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ABSTRACT 

We investigate the spatial properties of the large scale structure (LSS) of the Universe in the 
framework of coupled dark energy (cDE) cosmologies. Using the public halo catalogues from 
the Codecs simulations - the largest set of N-body experiments to date for such cosmolog- 
ical scenarios - we estimate the clustering and bias functions of cold dark matter (CDM) 
haloes, both in real- and redshift-space. Moreover, we investigate the effects of the dark en- 
ergy (DE) coupling on the geometric and dynamic redshift-space distortions, quantifying the 
difference with respect to the concordance ACDM model. At z ~ 0, the spatial properties of 
CDM haloes in cDE models appear very similar to the ACDM case, even if the cDE models 
are normalized at last scattering in order to be consistent with the latest Cosmic Microwave 
Background (CMB) data. At higher redshifts, we find that the DE coupling produces a sig- 
nificant scale-dependent suppression of the halo clustering and bias function. This effect, that 
strongly depends on the coupling strength, is not degenerate with <3% at scales r < 5 — 10 
h Mpc . Moreover, we find that the coupled DE strongly affects both the linear distortion 
parameter, p, and the pairwise peculiar velocity dispersion, On- Although the models con- 
sidered in this work are found to be all in agreement with presently available observational 
data, the next generation of galaxy surveys will be able to put strong constraints on the level of 
coupling between DE and CDM exploiting the shape of redshift-space clustering anisotropics. 
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1 INTRODUCTION 

The concordance ACDM cosmological model has strongly im- 
proved its reputation over the last two decades thanks to the 
overall agreement between its predictions and the e ver increas- 
ing wealth of available obse r vational data (see e . g. iRiess et alJ 
19981 : IPerlmutter et all 119991; IPercival et all |200lt ISpergel et al l 
20031: lAstier et alj|2006t iKomatsu et alj|2009| - ICrocce et alJl201lF 
Carnero et al l|). However, some possible tensions between the 



predictions of the ACDM scenario and astrophysical observations 
at small scales have also been claimed, like the abundance of satel- 
lite galaxies in CDM haloes, the low baryon fraction in galaxy clus- 
ters and the 'cusp-core' problem for the ha l o density profiles (see 
e.g. lAUen et al.ll2004 lLaRoque et al]|200fj: IVikhlinin et alj|200fj| : 
ISimon et al.ll2003l ; iBinnev & Evansll200ll : lNewman et alj|2009h . It 
is still unclear whether such discrepancies might be due to an incor- 
rect description of the baryonic phenomena at work, or if they are 
directly connected to the underlying cosmological scenario. In any 
case, it is interesting to investigate how the standard ACDM model 



could be modified in order to better match the above-mentioned 
observations. 

Furthermore, the fundamental nature of CDM particles and 
of the DE component - that together constitute more than 95% of 
the total energy of the Universe - is still unknown, which repre- 
sents one of the central problems of modern cosmology. From an 
observational point of view, the existence of a stress-energy com- 
ponent with negative pressure, responsible for the accelerated ex- 
pansion of the Universe, can be directly deduced from supernova 
measurements (see e.g. [Riess et all 1 1 99 8l; [Perlmutter et al.l IT999I : 
ISchmidt et ail 1 19981 ; lAstier et al J |2006t iKowalski et al.1 l2008h 7or 
by jointly combining the CMB power spectrum and low redshift 
observations, like measurements of the Hubble constant, of the 
galaxy power spectrum and of cluster counts. More recently, by 
combining the CMB lensing deflection power spectrum from the 
Atacama Cosmology Telescope with temperature and polarization 
power spectra from the Wilkinson Microwave Anisotropy Probe, it 
has been possible to break the geometric degeneracy and constrain 
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the da rk energy density usi ng measurements of the CMB radiation 
alone dSherwin et al]|201lh . 

Several theoretical efforts have been made in order to find a 
satisfactory explanation for the cosmic acceleration. Different DE 
models have been proposed in the literature, ranging from sim- 
ple scenarios with a constant equation of state to models con- 
sidering intera ctions between DE and other cosm ic fluids, as e.g. 
CDM (see e.g.lWetterichlll995l : lAmendoliJ 2000) or massive neu- 
trinos dAmendola et alj2008l) . Indeed, if the accelerated expansion 
is driven by a scalar field, there is no fundamental reason why DE 
and CDM should not interact with each other. Interestingly, it has 
been shown that such a direct coupling could help in alleviating 
so me of the small-scale problems of the ACDM cosmology (se e 
e.g. lBaldi et alfeOloLlBaldi & Pettoriiioboi lfclEee & BaldfeOl ll) . 

In the present paper, we investigate the spatial properties of 
LSS predicted by interacting DE cosmological models. In par- 
ticular, this work is focused on the clustering and bias proper- 
ties of CDM haloes and on their redshift-space clustering distor- 
tions (RSD) caused by the line-of-sight component of galaxy pe- 
culiar velocities. The latest galaxy surveys that have enabled mea- 
surements of RSD are the 2-degree Field Galaxy Redshift Su rvey 
dPeacock et alj200ll;lHawkins et al.l2003|:|Percival et alj20o4). the 
Sloan Digital Sky Su r vey dTegmark et alj2004IZehavi et alj2005l; 

Tegmark et in 12000I : Okumura et al. 2008; Cabre & Gaztanaga 

2009al lbh and the VIMOS-VLT Deep Survey dGuzzo et alj2008h 

RSD can be used to probe gravit y theories, a s they are directly 
related to the growth rate of LSS dLinderl 12008). Moreover, they 
can also be used in several other contexts, e.g. to robu stly constrain 
the va lue of the total mass of cosmological neutrinos dMarulli et al.l 
1201 lh or to investigat e the dynamical properties of the warm-hot 
intergalactic medium dUrsinoetalj201ll) . Since the interaction be- 
tween DE and CDM strongly affects the formation and evolution 
of LSS, a significant effect on RSD is also expected. To properly 
describe all the non-linear effects at work, we use a series of state- 
of-the-art N-body simulations for a wide range of different cDE 
scenarios - the CoDECS simulations dBaldill2011cl) - and discuss 
the impact that DE interactions can have on LSS with respect to 
the ACDM case. Both on-going and next generation galaxy sur- 
veys, such as VIPERS (Guzzo et al. 201 1 in preparation), BigBOSS 
dSchlegeletal .1120091) and Euclid dLaureiisll2009l ; iRefregier et all 
l2010l) . will greatly improve the precision in RSD measurements 
over an even grater redshift range, up to z = 2, allowing for direct 
constraints on the underlying cosmological model and on the grav- 
ity theory. 

The structure of the paper is as follows. In section ij2]we de- 
scribe the cDE models analyzed in this work, while in section Sj3]we 
introduce the exploited set of N-body experiments used to simulate 
the LSS in these cosmological frameworks. In section Sj4]we ana- 
lyze the CDM halo clustering and bias functions and the redshift- 
space clustering distortions. Finally, in section jj5] we draw our con- 
clusions. 



2 THE COUPLED DARK ENERGY MODELS 

Interacting DE scenarios have been introduced as a possible al- 
ternative to the standard ACDM cosmology dWetterictJ [l99l ; 
lAmendolj '2000) with the aim of addressing the fine-tuning prob- 
lems that characterize the cosmological constant. These alternative 
cosmologies include models where a DE scalar field <|) interacts 
with other cosmic fluids by exchanging energy-momentum during 
the evolution of the universe. Several possible forms of DE inter- 



action have been proposed in the literatu re, including e.g . mod- 
els of coupling to the CDM fluid (as e.g . [Amendola 2000, 120041 : 
iKovama et alj|2009l:lHonorez et alJbOldlBaldfol lot) or to mas- 
sive neutrinos dAmendola et alj|2008l)~ In the present work, we will 
focus on the former class, considering a DE-CDM interaction de- 
fined by the following set of background dynamic equations: 



dV 



p c + 3Hp c = 

p h + 3Hp h = 

p r + 4Hp r = 

3H 2 = 



3^W 



Pb 



(1) 

(2) 

(3) 
(4) 

(5) 



where an overdot represents a derivative with respect to the cosmic 
time t, H = d/a is the Hubble function, V(§) is the scalar field self- 
interaction potential, Mp\ =1/ \/%%G is the reduced Planck Mass, 
and the subscripts b,c ,r indicate baryons, CDM, and radiation, re- 
spectively. The strength of the interaction between the DE scalar 
field and CDM particles is defined by the coupling function T| c ((])), 
while the shape of the potential V((|)) dictates the dynamical evolu- 
tion of the DE field that in turn determines a time evolution of the 
CDM particle mass: 



d\nM r 



(6) 



dt V 3 

as one can obtain from Eq. 10. 

In this w ork, we will c onsider the set of coupled DE scenar- 
ios defined in iBaldil J201 lcl) . which includes two possible choices 
for the coupling function - namely a constant coupling r| f = T|o 
and an exponential coupling T| c ((])) = r|oexp [r|i<])]- and two possi- 
ble forms of the scalar potential, namely an exponential potential 
dLucchin & Matarresdl 19851 : IWettericbll 19881) : 



V(<S>)=Ae- ait > 
and a SUGRA potential dBrax & Martini 19991) : 



(7) 



(8) 



where for simplicity the field <|) has been expressed in units of the 
reduced Planck mass in Eqs. 17 1 8) . The parameters of the models 
investigated in the present work are summarized in Table Q] 



At the linear perturbations level, the time evolution of density 
fluctuations in the coupled CDM fluid and in the uncoupled bary- 
onic component (8 C j, = 5p c j,/p c .b, respectively) obeys the follow- 
ing system of equations: 



5 C 



-2H 



1-TV 



Hs/l 

l h = -2H8 h + 4%G [p b h h + p e 8, 



o c +47lG[pi5i + p c 5 c r c ] , (9) 

(10) 



where for simplicity the field dependence of the coupling function 
T| c (<])) has been omitted. In Eq. (9]l, the factor F c = 1 +4r|^((|))/3 
encodes the effect of the additional fifth-force mediated by the DE 
scalar field (]) for CDM perturbations. Additionally, the second term 
in the first square bracket at the right-hand-side of Eq. l[9} repre- 
sents an extra-friction term on CDM fluctuations arising as a con- 
sequence of momentum conservation (see e.g. lAmendol3 20041 : 
iPettorino & Baccigalupj|2008l : iBaldi et alJl2010l : lBaldill2011al for 
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Table 1. The list of cosmological models considered in the CoDECS project and their specific parameters. 
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-1.0 


0.809 


EXP001 


V($) =Ae 
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0.08 


0.05 





-0.997 


0.825 


EXP002 


V(4>) =Ae 


-ao 


0.08 


0.1 





-0.995 


0.875 


EXP003 


V((j>) =Ae 


-cut) 


0.08 


0.15 





-0.992 


0.967 


EXP008e3 


v\if)=Ae 


-CH|> 


0.08 


0.4 


3 


-0.982 


0.895 


SUGRA003 


V(4>)=Af 




2.15 


-0.15 





-0.901 


0.806 



a derivation of Eqs. dU5l9U0t and for a detailed discussion of the 
extra-friction and fifth-force corrections to the evolution of linear 
perturbations). 

The background and linear perturbations evolution of all the 
cosmo logies consid ered in this work h ave been discussed in full de- 
tail in lBaldil d2011bl) and lBaldil d201ld) , to which we refer the inter- 
ested reader for a more thorough description of the models. In par- 
ticular, the EXP008e3 model - based on an exponential potential 
and on an exponential coupling function - features a steep growth 
of the interaction strength at low redshifts and has been shown 
to significantly enhance the pairwise infall v elocity of collidin g 
galaxy clusters similar to the "Bullet Cluster" (Lee & Baldi 2oTlh . 
On the other hand, the SUGRA003 model based on a SUGRA po- 
tential and on a negative constant coupling is an example of the re- 
cently proposed "Bouncing" cDE scenario ( lBaldill2011bh , which is 
characterized by an inversion of the scalar field direction of motion 
at relatively recent epochs (z; nv ~ 6.8 for the specific model under 
investigation), giving rise to a "bounce" of the DE equation of state 
parameter wa on the cosmological constant barrier wa = — 1. As a 
consequence of such inversion, the quantity T| e <|) that appears both 
in Eqs. ® and (|9) changes sign in correspondence to the bounce 
such that the CDM particle mass decreases before the bounce and 
increases again afterwards, while the extra-friction term changes 
from a drag to a proper friction at Zim- Such peculiar dynami- 
cal behaviour has been shown to provide a possible explanation 
dBaldfcOllbh to the recent detection s of massive cluster s of galax- 
ies at high redshifts dJee et alj[2009l ; iRosati et al]|2009h . although 
the statistical significance of such detections as a challenge to the 
standard ACDM scenario is presently not ye t conclusive (see e.g. 
iMortonson et alj201ll ; IWaizmann etal.ll201ll) . 



3 THE N-BODY SIMULATIONS 

For our investigation we rely on the public halo catalogues of 
the CoDECS simulations (Baldi 2011c), the largest N-body simu- 
lations for cDE cosmologies to date, that include all the models 
described in Table Q] In particular, we make use of the L-CoDECS 
suite which consists of large scale, collisionless N-body runs fol- 
lowing the evolution of 1024 3 CDM particles and 1024 3 bary- 
onic particles in a cosmological box of 1 comoving Gpc//i on a 
side. T he simulations hav e been carried out with the modified ver- 
sion bv lBaldi et alj d2010h of the wi dely used parallel Tree-PM N- 
body code GADGET (Springel 2005), and have a mass resolution of 
m c (z = 0) = 5.84 x 10 10 M e /h and m b = 1.17 x 10 10 M Qj /h for 
CDM and baryons, respectively, and a force resolution of e g = 20 
kpc /h. 

The initial conditions for all the different cosmologies were 
generated by perturbing a ho mogeneous "glass" distribution 
dWhitdl 19941: iBaugh et aljfl995h based on the same linear power 



spectrum at zcmb> ar, d then rescaling the resulting displacements 
to the starting redshift of the simulations, zi = 99, with the spe- 
cific growth factor D + (z) of each model computed by numerically 
solving Eqs. l |9|10t . This procedure ensures that all the models are 
consistent with the same normalization of density perturbations 
at the last scattering surface, which has been assumed to corre- 
spond to the "WMAP 7 Only Maximum Likelihood" constraints of 
iKomatsu et al.1 d201 ll) . The cosmological parameters at z = as- 
sumed for all the models included in the CoDECS project are Hq = 
VOJkms-'Mpc- 1 , f2 C DM = 0.226, Q. DE = 0.729, o 8 = 0.809, 
Qj = 0.0451 and n s = 0.966. 

The halo catalogues have been obtained by identi fying groups 
of pa rticles by means of a Friends-of-Friends (FoF, iDavis et al.l 
Il985l) algorithm with linking length X = 0.2 x d, where d is the 
mean interparticle separation. The FoF algorithm was run over the 
CDM particles as primary tracers of the matter distribution, and 
then baryonic particles have been attached to the FoF group of their 
closest CDM neighbour. In order to identify CD M substructures 
we ha ve used the SUBFIND algorithm described in Sprin geTet al.l 
d200ll) . All the results presented in this paper have been obtained 
using mass-selected sub-halo catalogues, composed by the gravita- 
tionally bound substructures that SUBFIND identifies in each FoF 
halo. For all the simulations, we have restricted our analysis in the 
mass range M m j n < M < M max , where M mm = 2.5 • lO 12 M //z and 
Af max = 3.6- 10 15 , 1.1 • 10 15 ,4.9- 10 14 ,2.6- 10 14 , 1.8- 10 14 M e /7i at 
z = 0,0.55, 1, 1.6,2, respectively. For the analysis presented in this 
paper, we have used one realisation of each model. To investigate 
how reliable our results are, we have divided each simulation box 
in 27 independent sub-boxes and estimated the uncertainties in our 
predictions from the scatter between them. At the scales consid- 
ered, we found a very good agreement between the errors estimated 
in this way and the theoretical ones, both for the correlation func- 
tion and for the bias. 



4 LSS IN CDE MODELS 

In this section we show how the halo clustering, the bias function 
and the redshift-space distortion parameters (3 and O12 get modified, 
with respect to the standard ACDM case, when CDM and DE can 
interact with each other. 



4.1 The DM halo clustering in real space 

In cDE scenarios the formation and evolution of cosmic structures 
can differ significantly with respect to the ACDM case. To analyze 
the spatial properties of LSS as a function of the DE coupling, we 
measure the CDM halo two-point correlation function, %(r), de- 
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Figure 1. Upper panels: real-space t wo-point correlat ion function of CDM haloes for all the models of the CoDECS project. The black error bars represent 
the statistical noise as prescribed by IMo et alj 0992). Lower panels: percentage difference between cDE and ACDM predictions, A^ c qe = 100- (^ c de — 
^acdm)/5acdm- At £ = 0, the spatial properties of CDM haloes in cDE models are very similar to the ACDM case. However, at higher redshifts the DE 
coupling produces a significant scale-dependent suppression in the halo clustering. 



fined as: 



dP = n 2 [l+E,(r)}dV l dV 2 , 



(11) 



where dP is the probability of finding a halo pair with one of them 
in the volume dV\ and the other in the volume dV 2 , separated by 
a comovin g distance r. To mea sure the function %(r), we use the 
standard fLandv & Szalavl d 19931) estimator: 

HH{r) - 2RH(r)+RR(r) 



RR(r) 



(12) 



where HH(r),HR(r) and RR(r) are the fraction of halo-halo, halo- 
random and random-random pairs, with spatial separation in the 
range [r — dr/2, r + dr/2\. The random samples are three times 
larger than the halo ones. A computationally efficient linked-list 
algorithm has been used to count the number of pairs in the halo 
and random catalogues. 

Fig. Q] shows the real-space two-point correlation function 
of CDM haloes for all the models of the CoDECS project. The 
black err o r bars represent the statistical noise as prescribed by 
IMo et al] d 19921) . The lower panels show the percentage differ- 
ence between cDE and ACDM predictions, Ac^ c de = 100 • (cj C DE — 
^acdm)Aacdm- The z — spatial properties of CDM haloes in 
cDE models are very similar to the ACDM case. The EXP003 
model shows the largest deviations that are in any case quite small, 
between ~ 2% and ~ 6% for comoving separation < 50 hT x Mpc . 
The situation is different at higher redshifts, where the DE coupling 



produces a significant scale-dependent suppression in the halo clus- 
tering. This effect strongly depends on the coupling strength and 
on the redshift. At small scales, the deviations with respect to the 
ACDM predictions rise to ~ 20% at z = 1 and ~ 35% at z = 2 
for the extreme EXP003 model. The EXP008e3 and SUGRA003 
models predict a different redshift evolution in the clustering func- 
tion with respect to the models with a constant coupling. Therefore, 
measuring the redshift evolution of tj(r) can help to disentangle dif- 
ferent interacting DE cosmologies. 

4.2 The CDM halo biasing function 

The clustering suppression predicted by cDE models in mass- 
selected halo samples can be better understood by looking at the 
large-scale halo bias. Fig.[2]shows the mean apparent effective bias 
of CDM haloes defined as follows: 



< b(z) > = < 



Tihalo.cDE 



3DM.ACDM(o 8 ) 



(13) 



where ^ ha i 0iCDE is the cDE halo clustering and |dm,ACDM(o 8 ) is 
the ACDM DM clustering normalized to the CTs values of the 
CoDECS models. i; DM acdm(o 8 ) nas Deen obtained by Fourier 
trasforming the non-lin ear power spectrum extracted from CAMB 
Lewis & Bridle! |2002|) . which exploits the HALOFIT routine 
Smith et al. 2003). The bias is averaged in the range 10h~ l Mpc < 
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Figure 2. Mean apparent effective bias of CDM haloes, <b>=< (^halo.cDE/^DM.ACDM(o g ))°' 5 >■ ave raged in the range lOh ^ 1 Mpc < r <50hT l Mpc . Dashed 
and dotted black lines show the theoretical ACDM effective bias computed according to the relations of Sheth et al. (2001) and Tinker et al. ( 2010), respectively, 
normalized to the 08 values of the CoDECS models. The error bars represent the propagated statistical noise as prescribed bv lMo et all 119921) . while the grey 
shaded areas show a 10% error. The effect of DE coupling on the CDM halo large scale bias appears totally degenerate with o~g for all the cosmological models 
considered. 



r <50fc _1 Mpc. The er ror bars represen t the propagated statistical 
noise as prescribed by iMo et all dl992l) . By definition, < b(z) > 
represents the apparent bias that would be derived in a cDE uni- 
verse if a ACDM model was erroneously assumed to predict the 
DM clustering. 

Dashed and dotted blac k lines in F ig. [2] show the theoreti- 
cal AC DM effective bias of ISheth et all feoOll) and iTinker eTall 
(2010), respectively, normalized to the values of Cg of the CoDECS 
models and weighted with the halo mass function, n(M,z): 



b e s(z) = 



f$ m n(M,z)b{M,z)dM 



Su-n{M,z)dM 



(14) 



where M m ; n and M max have been defined in section ij3] The grey 
shaded bands show a 10% error, that approximately represents the 
uncertainty in the theoretical bias predictions presented in the liter- 
ature. 

The effect of DE coupling on the CDM halo large-scale bias 
appears totally degener ate with gg for all th e cosmological mod- 
els considered (see also ldemson et al1l201ll for a general discus- 
sion of degeneracies in cDE models). As already seen in Fig. Q~| 
the halo clustering properties at small scales can help in remov- 
ing this degeneracy. This can be better appreciated in Fig. [3] that 
shows CDM halo bias as a function of scale. The lower panels show 



the percentage difference between cDE and ACDM predictions, 
A£> c de = 100- (^cDE ~~ ^acdm ) /^acdm ■ Das hed and dot t ed lin es re- 
fer to the theoretic al ACDM predictions of Shet h et alj l l200ll) and 
ITinker et alj d2010l) , respectively, normalized to the values of Cg of 
the CoDECS models. It is clear that the scale-dependent suppres- 
sion of the halo bias at small separations, due to the DE coupling, 
breaks the Og degeneracy. Indeed, at scales r < 5 — 10 h Mpc , 
both the models with a constant coupling and the EXP008e3 one 
are not well reproduced by a standard ACDM model after having 
rescaled the matter power spectrum to the z = normalization of 
the cDE model. The SUGRA003 appears in better agreement with 
the ACDM predictions at z < 1, while at z > 2 the halo bias of this 
model is not degenerate with Og at all scales, although the effect is 
quite small (~ 3% at z = 2). 



4.3 Redshift-space clustering distortions 

4.3.1 Clustering shape in redshift-space 

The redshift of a galaxy does not correspond to a unique distance 
from the observer, as the line-of-sight component of a galaxy's pe- 
culiar velocity creates an additional Doppler shift. When the dis- 
tances are computed without correcting for this peculiar velocity 
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Figure 3. Upper panels: apparent effective bias of CDM halo as a function of scale, at three different redshifts z = 0, 1 , 2 from left to right. The error bars rep- 
resent the statistical noise )Mo et al.ll992l) . Lower panels: percentage difference between cDE and ACD M predictions , Afc c pE = 1 00 ■ (& c de — ^acdm ) / ^ACDM ■ 
Dashed and dotted lines show the theoretical ACDM predictions computed according to the relations of lSheth et alj feOOll) and Tinker et al. (2010), respec- 
tively, normalized to the z = values of Og of the CoDECS models. The scale-dependent suppression of the halo bias at small separations (r < 5 — 10 h Mpc ) 
due to the DE coupling breaks the 0"8 degeneracy. 



contribution, we say to be in redshift- space and we refer to the 
redshift-space spatial coordinates using the vector s (while f in- 
dicates real-space coordinates). A redshift-space galaxy map is dis- 
torted with respect to the real one. RSD provide crucial constraints 
on the build-up of LSS. On small scales (<1 h~ l Mpc ) the distor- 
tion is mainly caused by the random orbital motions of galaxies 
moving ins ide virialised s tructures, i.e. the well-known fingers of 
God effect Jjackson|[l972l) . On larger distance scales, the coherent 
bulk motion of virialising structures leads to an apparent excess in 
the clustering strength perpendicularly to the line-of-sight. 

Since both comoving coordinates and peculiar velocities are 
known in an N-body simulation, we can investigate the spatial 
properties of LSS both in real- and redshift-space. To obtain 
redshift-space halo maps, we locate a virtual observer at z = and 
put each simulation box at a comoving distance corresponding to 
its output redshifts. The redshift-space distance of each CDM halo 
is obtained through the following equation: 



L 



obs 



(15) 



where H(z' ohs ) is the Hubble rate, c is the speed of light and the 
observed redshift, z t, s , is given by: 



£obs 



Zc- 



V ll 



(16) 



Here z c is the cosmological redshift, due to the Hubble flow, while 
the second term of Eq. J 16b is caused by galaxy peculiar velocities 
(V|| is the component of the peculiar velocity parallel to the line-of- 
sight). 

Whether the cosmological parameters assumed in Eq. d 1 5 1 > are 
not the true ones, or the Universe is not described by the ACDM 
equations, like in the cDE scenarios, a geometric kind of distor- 
tion is added to the dynamic distortion caused by peculiar velocities 
dAlcock & Paczvnskill 19791) . However, as we will see in Sec. 14.3.21 
geometric distortions in cDE models are rather small relative to the 
dynamical ones and can therefore be neglected. 

Such a simple method to derive redshift-space coordinates is 
not as accurate as a full past light-cone reconstruction (see e.g. 
iMarulli et all2009l) , since it neglects the redshift evolution of halo 
properties inside each snapshot box. However, this method is accu- 
rate enough for the purpose of the present work. In a companion 
paper we aim to estimate the expected observational uncertainties 
in constraining cDE model parameters with next generation galaxy 
surveys, and for that specific goal realistic mock galaxy catalogues 
will be constructed. 

To extract cosmological constraints from RSD it is convenient 
to decompose the distances into the two components perpendicu- 
lar and parallel to the line-of-sight, r = (rj_ , r\\ ), When measured in 
real-space, the contour lines of £,(r^, rn ) are circles, as the popula- 
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tion of haloes/galaxies is isotropic when averaged on large scales. 
Instead, the redshift-space correlation function is distorted: at small 
scales, the fingers of God effect changes the shape of correlation in 
the direction parallel to the line-of-sight, while at large scales the 
bulk motion squashes the correlation perpendicularly to the line-of- 
sight. These effects can be clearly seen in Fig. [4] where the contour 
lines of the two-point redshift-space correlation function of CDM 
haloes are shown. We also plot in the same panels the undistorted 
real-space correlation function (dotted grey lines) for comparison. 

4.3.2 Modelling the dynamical redshift-space distortions 

Assuming the plane-parallel approximation, the linear power spec- 
trum of the matter density fluctuations, P(k), can be parameterized 
in redshift-space as follows: 

P(k) = (l + fr 2 ) 2 P lin (k), (17) 

where P\\ n (k) is the linear power spectrum in real-space, fj is the co- 
sine of the angle between k and the line-of-sight and (3 is the linear 
distortion parameter. Fourier transforming equation l !17b gives 

%{s,fi) = ^(^M+feWW + ^WW, (18) 

where the functions P; are the Legendre polynomials faamiltor] 
1 19921) . while the multipoles, *^n(s), can be written as follows 




6>W : 



(19) 



30 10 
s ± [Mpc/h] 

Figure 4. Contour lines of the two-point redshift-space correlation func- 
tion of CDM haloes, ^(i^.iy). The iso-curve plotted are ^(«x,'S||) = 
{0.05,0.08,0.13,0.2,0.3,0.5,1,3} and {0.08,0.15,0.25,0.5,1,3} in the 
left and right panels, respectively. Dotted grey lines show the undistorted 
correlation function in real-space. 



Us) 



4p 

3 



4p2 
7 



35 



%{r) + \%{r)-\t,{r) 



(20) 



(21) 



Here %{r) is the real-space correlation function and the barred cor- 
relation functions are defined as 



Jo 

■ fdrW)r' 4 
Jo 



(22) 



(23) 



To include in the model also the small non-linear scales, we con- 
volve it with the distr ibution function of random pairwise veloci- 
ties, /(v) (but see e.g. lScoccimarrdl2004lMatsubara1l2004h : 

S(*L,«||) = J~Jvf(vns ± ,s ll -v/H(z)/a(z)). (24) 

In this work, we adopt for /(v) the form 

1 / \/2> 



/exp(v) = 



<T 12 \/2 



exp 



0"12 



(25) 



where O12 is the dispersion in the pairwise peculiar velocities. 

To summarize, the model for the distorted correlation function 
in redshift-space (Eq. d24b ) depends on the correlation function in 
real space, ij(r), and on two free parameters: the linear distortion 
parameter, (3, and the pairwise peculiar velocities dispersion, 0"i2. 
It can be demonstrated that P is well parameterized by the following 
equation: 



f{a„ 



(26) 



where f(£l m ) = dlnD/dlna is the velocity growth rate, D is the 



linear density growth factor, b is the linear CDM halo bias, and 
fl m (z) is the matter density. The so-called gravitational growth in- 
dex, y, directly depends on the gravity theory, being y~ 0.545 for 
GR gravity. As a consequence, observational constraints on the pa- 
rameter (3 that can be obtained modelling the redshift-space cor- 
relation function directly translate to constraints on the 
gravity theory. 

Fig. [5] shows the ratio between redshift- and real-space cor- 
relation functions of the CDM haloes in all the CoDECS simula - 
tions. The error bars represent the statistical noise dMo et al.l 19921) . 
The black lines show the theoretical ACDM prediction given by the 
large-scale limit of Eq. d 1 8 b : 



1 - 



2p 

3 



5 ' 



(27) 



where P = — , and the effective bias b e g is given by Eq. J14b . 
Dashed and dotted black lines refer to the theoretical ACDM pre- 
dictions obta ined using the effective bias of She th et al. and 
lTinkeretal1(l2010h . respectively, normalized to the values of Os of 
the CoDECS models. The grey shaded area represents the propa- 
gated 10% theoretical bias error. 

As expected, the measured and model ratios %{s)/%(f) agree 
quite well in the ACDM model. The cosmic evolution of LSS in 
cDE models has the effect of significantly suppressing this ratio. 
However, as for the large-scale clustering and biasing function, 
RSD appears to be strongly degenerate with as at the scales and 
redshifts considered. Indeed, the mean ratio %{s)/%(f) of CDM 
haloes in the cDE simulation can be reproduced quite well by the 
expected ACDM prediction rescaled at the a% value of the cDE 
model. Only the SUGRA003 model appears not totally degenerate 
with ag at z = 0.5, but the effect is quite small. 

To extract the cosmological information from RSD both at 
small and large scales, we use Eq. d24t to model the redshift-space 
correlation function . where the real-space cj(Y) has been 
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Figure 5. Ratio between redshift- and real-space correlation function of the CDM haloes. Dashed and dotted black l ines show the theoretical ACDM predictions 
(Eq. )27t ) using the effective bias functions computed according to the relations of lSheth et all l20 0^)jind [ Tinker et al] l2010l) , respectively, normalized to 
the values of o"8 of the CoDECS models. The error bars represent the statistical noise jMoetalJll992l) . The DE coupling strongly affects the redshift- 
space clustering distortions. However such effect is strongly degenerate with 08 for the redshift considered. Only the SUGRA003 model appears not totally 
degenerate with Gg at z = 0.5, but the effect is small. 



measured using directly the comoving CDM halo model coordi- 
nates. The best-fit parameters, (3 and 0"i2, obtained with a standard 
% 2 minimization procedure are shown in the upper panels of Fig. [6] 
as a function of redshift. The errors on the best-fit parameter pre- 
dictions have been estimated dividing the simulation box in 27 sub- 
boxes, measuring (3 and O12 in each of them and rescaling their scat- 
ter by the square root o f the total volume of the simulation box (see 
e.g. iMarulliet al]|201lh . The strong effect of cDE on both these 
parameters can be appreciated in the lower panels that show the 
percentage difference between cDE and ACDM predictions. Both 
in the models with a constant coupling and in the EXP008e3 one, 
(3 and O12 are higher with respect to the ACDM case. This effect 
depends significantly on the coupling strenght and for the most ex- 
treme models considered raises to the level of ~ 20% for (3 and 
~ 40% for O12. Viceversa, the SUGRA003 model predicts a sup- 



pression of ~ 10% both in (3 and in a 12 at z ~ 0.5 — 1, while at 
higher redshifts the differences with respect to the ACDM model 
decrease. 

This result suggests that using an independent constraint on 
as, that for instance can be obtained from the clustering and bias 
function at small scales (see Sec. 14. l| and |4.2| l, the next generation 
of galaxy surveys will be able to put strong constraints on the cou- 
pling between DE and CDM exploiting the shape of redshift-space 
clustering anisotropies. 

4.3.3 Geometric distortions 

In cDE scenarios, the relation between redshift and comoving dis- 
tance is different from the one in the ACDM case. This means that 
if the ACDM equations are used to obtain galaxy maps from red- 



© 0000 RAS, MNRAS 000, 000-000 



Redshift-space distortions in cDE models 9 




250 - 



- e 



150 









ACDM — 




V EXP003 - 






1 \^ 





1 1.5 2 
redshift 



40 



- K 20 




- 



1 1.5 
redshift 





1 




0.8 








0.6 


1 






0.4 






<i 






0.2 









2 




1.5 


■t-°. 






1 


I 








b* 




< 


0.5 












— EXP001 




- Expooa 




- EXP003 







0.5 




EXP00Be3 
SUGRA003 




I 5 



- < - 



1 1.5 
redshift 



1 1.5 
redshift 



Figure 6. Upper panels: Best-fit parameters p(z) (left panel) and 012(2) 
(right panel) measured in all the CoDECS simulations. The error bars 
shown olny for ACDM represent the scatter in the measured p and Q\i 
obtained by dividing the simulation box in 27 sub-boxes, and rescaled 
by the square roo t of the total volume of the simulation box (see e.g. 
iMarulli et alj201 l|). Lower panels: percentage difference between cDE and 
ACDM predictions, Ap cDE = 100 ■ (P cDE - PacdmVPacdm (left panel), 
Ao cDE = 100 ■ (oi2,cDE - Oi2,acdm)/oi2,acdm (right panel). 



shift measurements in a Universe where DE and CDM interact, 
the correlation function appears distorted. To investigate the effect 
of suc h geometric distortions, also called Alcock-Paczynski distor- 
tions dAlcock & Paczvnski|[l979l) , we have repeated all the mea- 
surements described in the previous sections assuming the ACDM 
Hubble parameter H(z) in Eq. l |15t to measure redshift-space dis- 
tances also in the cDE simulations. Fig. [7] shows the percentage 
error on the best-fit parameters introduced in this way. In the three 
models with a constant coupling (EXP001, EXP002 and EXP003) 
these distortions slightly increase as a function of redshift, but only 
up to the level of ~ 1% for (3 and ~ 2% for O12 at z = 2, for the 
extreme EXP003 model. So, as expected, these distortions are quite 
small compared to the dynamical ones. Also in the EXP008e3 mod- 
els the geometric distortions appear clearly negligible with respect 
to the dynamical distortions. Only the SUGRA003 model shows 
large geometric distortions, at a level of ~ 5% for (3 and ~ 10% for 
O12, comparable to the level of dynamical distortions. This suggests 
that geometric and dynamic distortions on the redshift-space clus- 
tering of galaxies might be strongly degenerate for the "Bouncing" 
cDE scenario, and independent measurements as e.g. the redshift 
evolution of the halo mass function should be used to break the 
degeneracy. 



5 CONCLUSIONS 

In the present paper we have investigated the impact of a possible 
interaction between DE and CDM on the clustering and redshift- 
space distortions of CDM haloes, also in view of the large wealth of 
high-quality data expected from the next generation of large galaxy 
surveys. 

In cDE scenarios, in fact, the formation and evolution of LSS 



Figure 7. Percentage errors on the best-fit parameters p(z) (left panel) and 
012(2) (right panel) introduced if Eq. H5\ is used to measure redshift-space 
distances also in the cDE simulations. 



can be significantly different with respect to the ACDM case. The 
interaction induces a time evolution of the mass of CDM particles 
as well as a modification of the growth of structures determined 
by the combined effect of a long-range fifth-force mediated by the 
DE scalar degree of freedom and of an extra-friction arising as a 
consequence of momentum conservation. The combination of these 
effects can significantly alter the evolution of linear and nonlinear 
density perturbations as compared to the standard ACDM scenario, 
and provide a direct way to test this class of cosmological models. 
In particular, the clustering properties of CDM haloes from linear 
to highly nonlinear scales might show specific signatures allowing 
to distinguish a cDE universe from ACDM. 

In order to investigate the spatial properties of LSS in these 
cosmological scenarios it is therefore necessary to rely on a fully 
nonlinear treatment of structure formation in the context of cDE 
models. To this end, we made use of the public halo catalogues ex- 
tracted from the largest set of N-body experiments to date for such 
cosmologies, the CoDECS simulations. From such catalogues, we 
measured the clustering and bias properties of CDM haloes both 
in real- and redshift-space, as a function of the DE coupling and 
redshift. Moreover, we investigated the effects of the DE coupling 
on the geometric and dynamic redshift-space distortions, quantify- 
ing the difference with respect to the concordance ACDM model. 
We compared numerical predictions with theoretical expectations 
finding a good agreement in ACDM, as expected, and a strong de- 
generacy between the DE coupling and the matter power spectrum 
normalization, rjg. 

In particular, at z ~ the spatial properties of CDM haloes 
in cDE models appear very similar to the ACDM case, even if the 
CoDECS simulations are normalized to be consistent with the lat- 
est CMB data at last scattering. At higher redshifts instead, we find 
that the DE coupling produces a significant scale-dependent sup- 
pression in the halo clustering and bias function. Furthermore, the 
DE coupling strongly affects the redshift-space clustering distor- 
tions. We quantified this effect by modelling RSD in terms of the 
redshift-space correlation function and the distribution function of 
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random pairwise velocities and deriving the linear distortion pa- 
rameter, (3, and the pairwise peculiar velocity dispersion, O12, for 
all the CoDECS models. We find that for cDE models character- 
ized by an exponential self-interaction potential, both a constant 
coupling and an exponentially growing coupling determine higher 
values of p and O12 with respect to the ACDM case. On the con- 
trary, Bouncing cDE models based on a SUGRA self-interaction 
potential predict a suppression of P and O12 at z ~ 0.5 — 1, while at 
higher redshifts the differences with respect to the ACDM model 
decrease. 

We also analyzed the geometric distortions induced by the DE 
coupling, finding that these are always negligible with respect to 
dynamical ones for all the exponential potential models, while they 
might have a comparable effect as the latter for the Bouncing cDE 
scenario. 

Overall, we found a strong degeneracy of the DE coupling 
phenomenology with the amplitude of linear density perturbations 
defined by Os: the clustering, the bias function and the RSD in cDE 
scenarios can be well reproduced by a standard ACDM model after 
having rescaled the matter power spectrum to the Og normalization 
of each cDE model, for scales larger than ~ 5 — 10 /i _1 Mpc. 
However, we also find that the degeneracy is broken at smaller 
scales. Our results therefore suggest that it is possible to put strong 
constraints on the coupling between DE and CDM exploiting the 
shape of redshift- space clustering anisotropies and comparing the 
clustering and bias functions at small and large scales. 

To conclude, we performed an extended analysis of the im- 
pact of an interaction between DE and CDM on the clustering and 
redshift-space distortions properties of CDM haloes, using state- 
of-the-art N-body simulations for these cosmologies. Although the 
effect of the coupling is highly degenerate with at large scales, 
at smaller scales the degeneracy is broken, offering the chance that 
the next generation of wide and high-precision measurements of 
the galaxy distribution in real- and redshift-space might be able to 
tightly constrain or possibly detect a new interaction in the dark 
sector. 
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